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ABSTRACT 

In the neutron star, there is the relativistic degenerate electron gas which is the densest in 
the Universe. Due to electron-electron magnetic interaction, the differential susceptibility can 
equal or exceed 1, which results in the magnetic system of a neutron star being metastable or 
unstable states. The Fermi liquid of nucleons under the crust can be in the metastable state, 
while the crust is in the unstable state, which results in a phase transition occurring and the 
formation of domain structures in the later. The magnetizations of adjacent magnetic domains 
have opposite directions. The shearing stress acting on adjacent domains can result in a series of 
shifts or fractures in the crust. The crust releases the magnetic free energy, which corresponds to 
the bursts observed in magnetars. Simultaneously, a series of shifts or fractures which the deep 
crust closed to the Fermi liquid of nucleons undergoes can trigger the phase transition of the 
Fermi liquid of nucleons from the metastable state to a stable state. The free magnetic energy in 
the Fermi liquid of nucleons is released, which corresponds to the giant flares observed in some 
magnetars. 

Subject headings: star:neutron — magnetic fields — Pulsar: general 



1. Introduction 

In 1930, de Haas and van Alphen first ob- 
served that the magnetization M in normal met- 
als oscillates, under applied intense magnetic 
field and low temperature. The oscillatory func- 
tions are sinusoidal series with its fundamental 
frequency which can be described by the ex- 
tremal areas of the cross-section of the Fermi 
surface normal to the a pplied magnetic field 
(jLifshits fe Kosevich|[l95rl . Using the impulsive 
field method, Shoenberg found an unexpectedly 
high amplitude for the second harmonic and pro- 
posed magnetic interaction betw een the conduc- 
tion electrons ( Shoenbers] 1962 ). which was so 
called Shoenberg effect. He suggested that the 
magnetizing field is not applied field H but mag- 
netic induction B. When the differential magnetic 
susceptibility Xm — d(AnM) / dB exceeds 1 ( We 
use the Gaussian units in the paper) , the spatially 



uniform state of electron gas is thermodynamically 
unstable. Then, the electron gas rapidly evolves 
into a stable state, and has a spatially inhomo- 
geneous magnetic field wi th domains w hich are 
called as Condon domains (|Condonl ll966'). 

The magnetization and the de Haas- van Alphen 
oscillating effect for a relativistic degenerate 
electro n gas have been studied by many au- 
thors dVisvanathanl Il962t ICanuto fe ChIuVll968l 
Chudnovskvi Il98ll ). Previous studies about 
the magnetic susceptibility of neutron stars fo- 
cused on whether the observed field result s from 
a spontaneous magnetizat ion ( Lee et al. 1969t 
O'Connell fc Roussell Il97lh . This almost can- 
not occur because the neutron star is insuffi- 
ciently cool. The domain structure and the 
potential applications in neutron star's crusts 
have been postulated by Blandford and Wilkes 
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(iBlandford fc Hernquistlll982l:IWilkes fc Ingraham 
1989). The inhomogeneous field resulting from do- 
main formation in the crust of a normal neutron 
star is negligible and can not produce observable 
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effects. However, for magnetars with magnetic 
fields stronger than normal neutron stars, the in- 
homogeneous field is significant so that magnetars 
can produce super-Eddington X-ray outbursts. 

Magnetars including Anomalous X-ray Pulsars 
(AXPs) and Soft Gamma Repeaters (SGRs) are 
characterized by their inferred dipolar magnetic 
field strength ranging from 5.9 x 10 13 G to 1.8 x 
10 15 G and their long spin periods ranging from 
5.2 s to 11.8 s. Surprisedly, some magnetars have 
undergone giant flares (GFs) in which the energy 
up to 10 46 ergs is released in a fraction of a sec- 
ond via 7- ray emissions. Up to now, three GFs 
ha ve been observed from respective ly SGR 0526 



66(lCline et alJll982ft . SG R 1900+14(lMazets et al 
1999J) and SGR 1806-20 ^Palmer et al.ll2005h . 



Their X-ray and 7-ray luminosities during 
persistent and burst phase are too large to be 
powered by their kinetic energy. Where is the 
magnetic energy released stored prior to the 
GFs? Is it i n the magnetospheres or in the 
neutron star? Duncan fe Thompson! (1992) and 
Thompson fc Duncan! (|2001 ) proposed a popular 
model for magnetars, and they considered that 
the released magnetic energy is stored in the neu- 
tron star. The controversy for the popular model 
is about it s triggering mechanism for h igh energy 
radiation. Duncan fe Thompson! ( 1992 ) suggested 
that a helical distortion of the magnetic field in 
the core induces a large-scale fracture in the crust, 
and a twisting deformation of the magnetic field in 
the crust and the magnetospheres. A GF may in- 
volve a large disturbance which probably is driven 
by a rearrange ment of magnetic field in th e deep 
crust and core (|Thompson fc Duncanlll995h . while 
the persistent emission can be explained by the 
very slow transport of t he field from the core to 
the c r ust by Hall drift (iGoldreich fc Reisenegger 
1992( 1. iKondratvevI j2002j) first postulated the ex- 
istence of magnetic domains in the magnetar's 
crust because of the inhomogeneous crust struc- 
ture. He argued that the burst activity of SGRs 
could originate from magnetic avalanches. How- 
ever, because of the domain forming mechanism 
resulting from the interaction between nucleus 
spins, the required magnetic field is in the range 
10 16 ~ 10 17 G. Whether the magnetic field in 
the crust of neutron star can be such strong is 
still controversial. Having considered the GF's 
emission energy (10 42 — 10 46 ergs) and it's mean 



waiti n g time dMazets et ah 1979t Hurley et al 
19991 : iTerasawa et all l2005h . IStella et all (|2005l l 
estimated that the internal field strength of mag- 
netar at birth time can reach up to B > 10 15 7 G. 
Later, having numerically calculated the magnetic 
properties of magnetar-matter such as the magne- 
tization and the susceptibility of elec tron, proton 
and neutron, ISuh fc Mathews! (|2010f l suggested a 
magnetic domain model to correlate smoothly be- 
tween the statistics of star-quakes and magnetic 
avalanches in magnetar crusts. 

The phase transition to the domain phase in 
magnetars is different from the ferromagnetic 
phase transition. The magnetic ordering in iron 
comes from the exchange interaction of bound 
electron spins at a sufficient low temperature, and 
it is not prerequisite to have an external magnetic 
field. While the magnetic ordering in magne- 
tar arises from the interaction between the or- 
bital magnetic moments of free electrons under 
high-quantizing field. Therefore, given a nearly 
uniform distribution of electrons, the phase tran- 
sitions in magnetars are similar to those in beryl- 
lium whe re magnetized m atter is the conductive 
electrons! Shoenberdll962l ). The relativistic Fermi 
sea of electrons is only slightly perturbed by the 
Coulomb forces of nuclei in the crust of magnetar 
and does not efficiently screen nuclear charges. 
The dominant contribution to differential suscep- 
tibility in a magnetar comes from electrons, while 
the contributions from nucleons are negligible (the 
magnetic moments of nucleons are three orders 
lower than electron orbital moments). 

In this paper we analytically calculate the dif- 
ferential susceptibility of relativistic degenerate 
electron gas, and find it is an oscillated function 
of magnetic field which is called as the de Haas 
and van Alphen oscillation. In §2, we discuss the 
magnetic phase transition. The differential sus- 
ceptibility is calculated in §3. In §4, we use our 
model to explain the observed evidences in mag- 
netars. Summary is in §5. 

2. The magnetic phase transition 

The magnetization of degenerate electron gas 
at high-quantizing magnetic field and low temper- 
ature [KT <SC hco c ) should exhibit a nonlinear de 
Haas- van Alphen effect. Without the magnetic in- 
teraction the oscillated magnetization are assumed 
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to be p eriodical in magnetic field H ra ther than 
in l/H (IPiDDardlll98nt lShoenbere l ll984l ): 



M = M sxn(H/H Q ), 



(1) 



where Hq hardly varies over one cycle of oscilla- 
tion. However, if we consider the feedback con- 
tribution from 4ttM, the magnetization depends 
on not only the quantizing magnetic field but also 
the cooperating ordering of the magnetic moments 
themselves. The magnetization actually should be 
a function of magnetic induction, M — M(B). 
Then, the magnetic induction is given by B = 
H + AttM(B). Replacing H by B is called as the 
Shoenberg or H-B effect. The oscillated sinusoidal 
function of the magnetization on magnetic field in 
Eq.(l) is given by 

M = Mo sm(B/H a ) = M sm[(H + AttM)/H ], 

(2) 

where Mq and Hq are constants in a period. 

When AttMq/Hq > 1, as Fig.l shows, M be- 
comes a three- values function of H, and the differ- 
ential magnetic susceptibility Xm = d(4irM) / dB 
can exceed unity. Having considered the case of 
a thin rod oriented along the direction of mag- 
netization, iPippardl ((1963) showed that the multi 
valued function in Fig.l is not physical. In the 
region of the curve between L and L' where the 
slope is lower than unity, the magnetized states 
are unstable and can never really exist. For any 
weak perturbation 5H near Hq, the perturbation 
of magnetic induction SB is given by 



SB 



5H 

1 - Xn 



(3) 



Obviously, there is a singular point when Xm = 1 
in Eq.(3), where the magnetized state is unstable 
and the first-order phase transition should occur. 
Then, the magnetic system should be in a stable 
state in which magnetization is inhomogeneous 
and magnetic domains form. The magnetizations 
in the adjacent domains have opposite directions. 
The stable magnetized states are represented by 
the dashed line between N and N' in Fig.l. But if 
there is a surface energy at the boundary of two 
different magnetization it may become metastable, 
which is similar to supe rheating or supercooling 
in gas- liquid transition ( Reich] 19981) . The solid 
line between L and N or L' and N' in Fig.l rep- 
resents the metastable state. It is not difficult 



to achieve the condition of the above magnetic 
phase transition in the terrestrial lab. The exper- 
iments of Condon indicated that the magnetized 
phase transitio n can take place in some metals 
( Condon Il966t ). However, what is the situation 
in compact objects such as neutron star? In the 
following we calculate the differential susceptibil- 
ity of relativistic degenerate electron gas and show 
the observable effects of magnetic phase transition 
in neutron star. 



3. The differential susceptibility of rela- 
tivistic degenerate electron gas 

The assembly of electrons in neutron star 
under a strong magnetic field is degenerated 
and relativistic. The ene rgy eigenvalues are 
( Johnson fc Lippmann 1949f ) 



E = [c 2 pl + /i 2 + /ie c (2n + s + l)] 1 ' 2 



(4) 



where fi = m c c 2 and e c = heB/cm c = hio c are 
the rest energy and cyclotron energy of electron, 
respectively. Here, p z is the momentum compo- 
nent along field direction which is the z-direction, 
n = 0, 1, ... and s = ±1 are the Landau and spin 
quantum number, respectively. The density of 
states per unit volume for the energy level is given 
by (pz/h)(eB/hc) = g(p z /h) where g = eB/hc is 
the density of states per unit area for a Landau 
energy level. As a consequence of quantization, 
the spherical Fermi surface of the free electron is 
replaced by a set of circles located on a spherical 
surface with a common axis along the B direction, 
which is shown in Fig. 2. 

The susceptibility of the electron assembly can 
be obtained by finding its grand potential which 
depends on the density of states. It is a function 
of energy and magnetic induction. The density of 
states per unit volume is given by 

n,s 

(5) 

where e = E — \l is the kinetic energy of an elec- 
tron. The grand potential per unit volume of the 
assembly is given by 



J = -/3/ln[l + e' 9 ^- e )]dZ(e,S) 
= -f °°Z(e 1 B)f(e)de, 



(6) 
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Fig. 1. — Magnetic field H vs. magnetization 
AttM 




Fig. 2. — The free-electron Fermi surface in the 
presence of a magnetic field along the p 2 -axis. 



where (3 = (fcT) -1 , if) is the chemical potential, 
/(e) is the Fermi-Dirac distribution function. Af- 
ter summing on spin quantum number, Eq.(5) re- 
duces to 



Z = ^3/2^/2 



where 



b = 



lb] 

2j^(b-n)W 



2e/j, 



2e c fj, 



(7) 



(8) 



and [b] is the integer of b ([b] < b). For a neutron 
star, the electron system is almost completely de- 
generate (ip ^> and the Fermi-Dirac distri- 
bution function is almost a step function except 
the region near e = if). The magnetic moment 
depends on the first-order derivative of grand po- 
tential with B and the dominant contributions 
come from the integral of Eq.(6) at the region 
near e = if). Hence, we define that b m — b\ t=1 p. 
If the chemical potential of neutron star is ~ 10 
MeV and the magnetic field is the quantum field 
Bq = 4.414 x 10 13 G (the cyclotron energy of elec- 
tron equals to its rest energy), then [b m ] is about 
200. Based on Eq.(7), we can see that the first- 
order derivative of grand potential with B has the 
singularity when b m = [b m ]. Therefore, the sus- 
ceptibility of rclativistic degenerate electron gas 
can be equal to or even larger than 1. 

For revealing the oscillating effect of the grand 
potential, we calculate the su m of Eq.(7) w ith the 
Poisson summation formula (|Dingldll952j ): 



oo oo „[(,] 

-F(0) + J2F(n)= ]T / F(x)e i2 ™dx. 

n—1 r— — oo 

(9) 

Because [b] is much larger than 1, the density of 
states in Eq.(5) is approximately by 



_jL 3/2r2L3/2 , y-oc 
1 v^oo 



{-iy' z {2v-l)\B 



(■/ + 2)/2 



2^2- 



2,4 V5(46)"(i/-1)!(j/+2)! 
cos(27rrb-37r/4) 1 

^3/2 \ J 



(10) 

where B m denotes the Bernoulli numbers. The 
first two terms on the right-hand-side are the non- 
oscillating parts while the third term is the oscil- 
lating part. 



The non-oscillating grand potential per unit 
volume can be obtained from Eq.(10). As the first- 
order approximation, we can keep only the first 
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term in the summation and obtain 
J= ~ 3^kr<^ /2 Jo°V + ^) 3/2 f(e)de+ 

6v^ ^ JO ( e 2 +2eM )i/2 ■ 

(11) 

In general, -0 3> [i for neutron star. Using the 
standard methods of evaluating integrals for the 
Fermi-Dirac distribution, we can evaluate the inte- 
grals in Eq.(ll). The non-oscillating susceptibility 
is given by 



and 



Xn 



= 4-7T 



d 2 J 



aw 

4V2^ ff 3/2 [ln M 



7^ /kT\2] 



(12) 



In our work, the chemical potential and the mag- 
netic field in the deep crust of neutron star is 
about 10 MeV and 10 12 G, the non-oscillating 
susceptibility can be ~ 5.6 x 10~ 3 . 

Using a approximate formula for the Fermi- 
Dirac distribution 



o Jo 



and the Fresnel Integrals 
J^cosCf^d* 



(13) 



x 
i 

2 



[x < 1) 

^sinCfjc 3 ) (z»l), 
(14) 

we can obtain the oscillating term of the grand 
potential per unit volume as 



J 



V2 3/2 V^oo 



iff 



y-oo 1 r/1 



_ ^(fcT) 2 (a r 0) 2 ] sin(a^ 2 



where 



An 



«(-) a (^-) 



V\ 3 ^ J> \-3/2 



(17) 



2^3 > 
3 V ' 



(18) 



Here, a in Eq.(17) is the fine structure constant. 
The result of Eq.(1 6) is similar with the result in 
IVisvanathanl ( 19621 ). It indicates that the sus- 
ceptibility oscillates with 1/B. We replace B 
by B' = B - B . When B -> B Q we have 
1/B « l/B (l - B'/Bq). Note that the first term 
is a constant, thus the oscillation are the period 
functions of B' . The differences of magnetic in- 
duction in an oscillating period are the periods of 
B', which are presented by 



SB 



2hec 2 B 2 
ip 2 r 



(r = l,2,-..). (19) 



jj) cos(a r /z 2 + 1 7r)+ 



(15) 

where a r = 7rr/(e c /i). On calculating the differ- 
ential susceptibility which is the second derivative 
of Eq.(15), we only keep the most rapidly vary- 
ing terms, that is, only differentiate cosines and 
sinusoidal. Because ip 3> /i for neutron star, the 
last term in Eq.(15) is dominant. The oscillating 
susceptibility is approximately given by 

^ ra= yl g( r -V 2 _A 1 (^)V/ 2 )co S (a^ 2 -^) ) 

r— 1 

(16) 



The coefficients of the cosine functions in Eq.(16) 
include two terms. The first term is the result of 
complete degeneracy, while the second term comes 
from the thermal fluctuation and is proportional 
to the second power of temperature. For a typi- 
cal neutron star, chemical potential, magnetic field 
and temperature are the order of 10 MeV, 10 12 G, 
and 10 6 K, respectively. The second term can be 
ignored because of KT <C tvjj c <C ipo- The oscil- 
lating susceptibility, Xm> is mainly determined by 
A . Fig. 3 shows the relation between logA and 
B/Bq. Obviously, the differential susceptibility 
can equal or exceed 1 when B <~ 15Bq. In this 
work, we assume that magnetars are similar to 
normal neutron stars except magnetic field: for a 
normal neutron star, B = 10 12 G while B = 15Bq 
for magnetars. For a stronger magnetic field, the 
condition fouj c <C ipo or [&m] 3> 1 cannot be sat- 
isfied, the above approximate methods cannot be 
used. We will discuss them in the next paper. 

4. The observed evidences in magnetars 

As mentioned in the last section, the suscep- 
tibility of relativistic degenerate electron gas can 
be equal to or exceed 1. The electron gas under 
strong magnetic field in neutron star may be in 
unstable state, and the first-order phase transition 
should occur. Finally, the electron gas should be 
a stable state and has the Condon's domain struc- 
ture. The magnetizations in the adjacent domains 
have opposite direction. 

Based on Eq.(19), the relative variation of mag- 
netic induction between adjacent magnetic do- 
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mams is 



SB ^ 
B Bq 



(20) 



For a normal neutron star(£?o ~ 10 12 G), SB /Bo ~ 
10~ 4 , it is very difficult to observe the oscillatory 
effects. However, for magnetars (Bo ~ WBq), 
SB I Bo ~ 10~ 2 , which is large enough to produce 
the bursts in SGRs or AXPs. 

The adjacent magnetic domains have different 
densities of electrons. The different magnetiza- 
tions in the adjacent magnetic domains mean that 
the phase difference is w. According to Eq.(16), we 
have S(a r i/j 2 ) = it. For simplicity we approximate 
the chemical potential tp with zero temperature 
free-field Fermi energy ip — chfin/Sn) 1 / 3 . The 
difference of electron density between the adjacent 
magnetic domains is given by 

Sn , e r . „ , Bn . 

^ W' (21) 

where e = ipo/p and we have taken fixed B = Bo- 
For magnetar, tpo ~ 10 MeV, Bo ~ 15£?q, and the 
rest energy of electron \x ~ 0.5 MeV. Obviously, 
here, e is scaled to 20 and B to 15Bq, respec- 
tively. 

The electrons in a neutron star coexist with 
other components such as nuclei in the crust or 
Fermi liquid of protons and neutrons closed to the 
deep crust. Due to the Coulomb interactions be- 
tween electrons and nuclei or protons, the mat- 
ter of neutron star may not be homogeneous once 
the domain structure appears in the crust. How- 
ever, it can not occur in the Fermi liquid. The 
mechanism of forming electron domain is the inter- 
actions between orbital magnetic moments which 
are counterbalanced by the degenerate pressure of 
electrons P c . Because the Coulomb interaction of 
the crystal lattices in the crust is the same order 
of magnitude with the interaction between orbital 
magnetic moments, the electron domain structure 
can lead to the forming of the nuclei domain struc- 
ture. However, in Fermi liquid of nucleons, the 
neutron degenerate pressure P n >• P c — Y e P n , 
where Y e is the fraction of electrons per neutron 
and it is only several tenths. The interaction be- 
tween orbital magnetic moments can not coun- 
terbalance against the neutron degenerate pres- 
sure P n . Therefore, the magnetic domain struc- 
ture cannot form in Fermi liquid of nucleons. The 
magnetic system is in a metastable state, which is 



similar to the supercooling or superheating states 
in the first-order phase transition. 

In the crust of magnetar we can roughly calcu- 
late the size of the domain structure. With the 
density increasing, the quantum number of Lan- 
dau energy also increases. Using the gravitational 
potential energy of a nucleon and the Fermi en- 
ergy per electron, we can estimate the heig ht of 
magnetic 



electron, we can estimate tne hcigJ 
domains by (|Suh fc MathewdbOlOh 



B 



V 0.36 A 20 ; V 15B Q ' 



(22) 

where g is the surface gravity and Y e is scaled to 
0.36 when neutrons drip out from nuclei. If we as- 
sume that the volume change of the electron gas 
during forming domains is only perpendicular to 
the radial direction of neutron star and the chang- 
ing magnitude (or the thickness of domain walls) is 
~ a cyclonic radius of electron, based on Eq.(21), 
the width of magnetic domain can be estimated 

f-™-'®'^ 2 cm ' (23) 



The Maxwell shearing stress between the adja- 
cent domains can deform the crust and result a 
strain, 9. It is given by 



B 6B 

47T 



91% 
4tt : 



(24) 



where v is the shear modulus of the crust and 
B v = \J Airv. Clo sed to the bottom o f the crust, 
B„-6x 10 15 G (|Bavm fc Pineslll97ll ). For mag- 
netars, Bo ~ 15-Bq, according to E qs.(20) and 
(24) , we have the strain 9 ~ 10 



Ruderman 



(|1991I ) calculated that the maximum strain of the 
crust, 6* max , is ~ 10~ 2 — 10~ 4 . Our result is in its 
range. 

When the domain structures appear in the crust 
of magnetar, the shearing stress acting on adja- 
cent domains can produce a relative shift or the 
fracture between them. This sudden shift or frac- 
ture can propagate with Alfven velocity along the 
domain layers, which results in a series of shifts 
or fractures in magnetic domains. The time scale 
of shifts propagating along a domain layer is esti- 
mated by 



I 

Va 



0.1( 



B P . . I 



15S Q 



)-*(_£_)( J_) 



(25) 
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wher e Va and p are the Alfven velocity ((Duncan 
120041 ) and the matter density at the deep crust, 
respectively, I is the length of domain layer and 
i?* (~ 10km) is the radius of magnetar. The time 
scale is agree with the observed dur ation of SGRs 
bursts ((Thompson fc Duncan! Il995l ) . The avail- 
able magnetic free energy is then 



5z_ 
10 3 " 



ergs, 



(26) 



where we scaled Sz to 10 3 cm (See Eq. (22)) and 
the SB to 10 13 G (See Eq.(20)). This energy also 
agree with the observation of bursts. 

The Fermi liquid of nucleons in magnetar may 
be in metastable state. A series of shifts which 
the deep crust closed to the Fermi liquid of nucle- 
ons undergoes can trigger the phase transition of 
the Fermi liquid of nucleons from the metastable 
state to a stable state. The free magnetic energy 
in the Fermi liquid of nucleons is released and it 
is far greater than that in the crust because the 
magnetic induction and the thickness of the Fermi 
liquid of nucleons is larger than the crust. Let the 
magnetic induction of Fermi liquid of nucleons be 
at the same order of magnitude with that in deep 
crust, and the thickness of Fermi liquid of nucle- 
ons be approximated by the radius of magnetar. 
We estimate that the energy released is 



1° 44 2 "^ 



(27) 



This energy is agree with those released in GFs of 
some magnetars. 



5. Summary 

We discussed the magnetization effects of rel- 
ativistic degenerate electron gas in neutron star. 
Having considered the magnetic interaction be- 
tween electrons, we found that the magnetic sys- 
tems may be unstable or metastable when the dif- 
ferential susceptibility equals or exceeds 1. Un- 
der an ultra-strong magnetic field, the magnetic 
domain structures in magnetar can appear in the 
solid crust, while Fermi liquid of nucleons may be 
in metastable state. The shearing stress acting on 
adjacent domains can result in a series of shifts or 
fractures in the crust. The crust releases the mag- 
netic free energy, which corresponds to the bursts 



observed in magnetars. Simultaneously, a series of 
shifts or fractures which the deep crust closed to 
the Fermi liquid of nucleons undergoes can trigger 
the phase transition of the Fermi liquid of nucleons 
from the metastable state to a stable state. The 
free magnetic energy in the Fermi liquid of nu- 
cleons is released, which corresponds to the GFs 
observed in some magnetars. 
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